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ABSTRACT 

We give a rigorous derivation of a theorem showing that charged particles in an arbitrary electromag¬ 
netic field with at least one ignorable spatial coordinate remain forever tied to a given magnetic-field 
line. Such a situation contrasts the significant motions normal to the magnetic field that are expected in 
most real three-dimensional systems. It is pointed out that, while the significance of the theorem has not 
been widely appreciated, it has important consequences for a number of problems and is of particular 
relevance for the acceleration of cosmic rays by shocks. 

Subject headings: Cosmic rays: general — magnetic fields — diffusion — particle acceleration — shock 
waves 


1. INTRODUCTION 

The problem of particle motion in electromagnetic fields 
is fundamental and pervades both laboratory and astro- 
physical plasma physics. In spite of the conceptual sim¬ 
plicity of the problem, theoretical analysis and numerical 
modeling of the particle motions can be difficult. Often, 
in order either to simplify the analytical work or to reduce 
the demand on computer resources, the system of inter¬ 
est is considered in only one or two spatial dimensions, 
with the hope that no essential physics will be lost. It is 
readily shown that this approximation does have a num¬ 
ber of important general consequences which can result 
in the loss of essential physics, and it must therefore be 
used with extreme caution. On the other hand, some real 
physical systems have a symmetry such that one (or even 
two) coordinates are ignorable. In such cases this imposes 
a real constraint on the possible motions of charged parti¬ 
cles. While we will discuss an example of this type of real 
symmetry, the primary emphasis of this paper will be to 
address the profound effect that this reduction of dimen¬ 
sionality has on theoretical stud ies of particle accele ration. 

Jokipii, Kota and Giacalone (Jokipii et al. 1993), here¬ 
inafter called JKG, presented a general theorem regard- 
ing spatial constraints to charged-particle motion in an 
electromagnetic field which has at least one ignorable spa¬ 
tial coordinate, and discussed explicitly the application of 
this theorem to hybrid simulations of collisionless shock 
waves. In the context of hybrid plasma simulations the re¬ 
sult has been stated in previous papers (Thomas & Brecht 
1988, Thomas & Winske 1991), so it has been known for 


a while by some workers in the field. In auroral and ra¬ 
diation belt physics it has been known and applied for 
much longer ( [Stormer 1955| ). The essence of this theorem 
is that in such systems a charged particle is effectively 
forever tied, in a sense which we will define later, to the 
same magnetic line of force, except for motion along the 
ignorable coordinate, or for the case of a magnetic field 
entirely aligned in the ignorable direction. Cross-field mo- 
tion is expected in many r eal, three-dimensional systems 
(Giacalone & Jokipii 1994). This constraint on particle 
motion represents a critical loss of physics in many sit¬ 
uations in which the ignorability of the coordinate is an 
artifact of the analysis rather than a property of the real 
system. Since JKG presented the theorem in a heuristic 
manner, and, because it appears that its implication has 
not been always understood or believed, in the present pa¬ 
per we present a rigorous derivation of the theorem and 
include an expanded discussion of its implications. 

The most salient application of the theorem concerns 
one- and two-dimensional hybrid and full plasma simula¬ 
tions, which have been employed by numerous researchers 
to study particle acceleration in collisionless shocks. A 
general conclusion of these investigations is that while the 
simulations can generate substantial non-thermal popu¬ 
lations of particles in quasi-parallel shock systems, none 
have been able to demonstrate any significant acceler¬ 
ation of a stochastic nature of ions and electrons in 
quasi-perpendicular shock environments. Hybrid simula¬ 
tions that model the quasi-parallel portion of the Earth’s 
bow shock (Trattner & Scholer 1991, Scholer, Trattner & 
Kucharek 1992) produced nice comparisons with AMPTE 
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suprathermal ion observations, following hard on the heels 
of the successful modeling of this data by Ellison, Mobius 
and Paschmann (1990) using a Monte Carlo technique to 
describe particle convection and diffusion. In contrast, 
highly oblique and quasi-perpendicular hybrid simulation 
models of the acceleration of interstellar pick up ions at 
the heliospheric termination shock (e.g. Liewer, Goldstein 
and Omidi 1993, Kucharek and Scholer 1995, Liewer, Rath 
and Goldstein 1995) fail to create any measurable non- 
thermal ions for field obliquities above around 50°. For 
full plasma simulations, which have had only very limited 
application to shock acceleration problems, the situation is 
similar. It has been concluded by such studies ( |Gallant~e^ 
al. lf)|92|) that perpendicular relativistic electron-positron 
shocks, such as would occur in pulsar wind termination 
shocks, are incapable of accelerating the ambient parti¬ 
cles to non-thermal energies. This situation is modified 
slightly by the inclusion of a proton component (Hoshino 
et al. 1992) that spawns magnetosonic waves that can 
provide limited (second order stochastic) acceleration of 
positrons only. We suggest that, in light of the following 
discussion, such results may be viewed as the consequence 
of the reduced dimensionality and the concomitant loss of 
cross field diffusion. 


2. THE THEOREM, ITS SPECIFICATION AND MEANING 

The theorem states that for a system of electromagnetic 
fields that have at least one ignorable coordinate the corre¬ 
sponding component of canonical momentum is conserved. 
The conservation of this component can have serious im¬ 
plications for the allowed motions of particles in such a 
field system e.g. binding a particle to a given field line 
equivalence class (FLEC) which is defined to be the class 
of all magnetic field lines that differ only by translation 
along the ignorable coordinate (or a rotation if the ignor¬ 
able coordinate is an angle variable). 

In the Appendix we show that if the electromagnetic 
fields do not depend on at least one coordinate, say X 2 , 
in some orthogonal (not necessarily Cartesian) coordinate 
system, then a particular gauge can be found in which the 
vector and scalar potentials, A and <p, also do not depend 
on this coordinate. Furthermore, insofar as the compo¬ 
nent A .2 is concerned, this gauge is unique up to an overall 
constant in space and time. This implies that the Hamilto¬ 
nian of a particle in such a field will have X 2 as an ignorable 
coordinate. 

We will also show in Section [2.11 that the scaled com¬ 
ponent of the vector potential, ^ 2^2 (a scale factor, hi, 
converts a generalized orthogonal coordinate, Xi , into con¬ 
ventional length units such as centimeters) in the ignor¬ 
able direction is a constant along a magnetic field line (as 
well as along the X 2 coordinate) and thus may be used as 
a label for a given FLEC where a field line equivalence 
class is defined as all field lines that differ only in their 
X 2 coordinate i.e. a surface generated by moving a field 
line through all values of the coordinate X 2 ■ We will also 
show that a physically meaningful magnetic field line mo¬ 
tion may be defined. This motion reduces to the usual 
convection of field lines for ideal hydromagnetic systems 
and can not be due to a simple gauge relabeling since no 
remaining gauge transformations that preserve the sym¬ 
metry can change the value of / 12 A 2 . We will also show in 


Section 
order of 


2.2 that the physical scale of this binding is of the 


the harmonic mean of the particle’s gyroradius in 


this field region. 


2.1. The equivalence class and its motion 

Since we have shown that with our choice of gauge noth¬ 
ing depends on the X 2 coordinate we may establish an 
equivalence class of all field lines that differ only in their 
X 2 coordinate. Further, from the definition of B, and the 
fact that BxB = Bx(\'xA) = 0, we may take the 2 
component of B x B and obtain 


Bi 


1 d(h 2 A 2 ) 

hi dxi 


+ B 3 


1 <9(112^2) 
h 3 dx 3 


B-V{h 2 A 2 ) = 0 (2-1) 


and / 12 A 2 is constant along a magnetic field line. This 
means that we may identify a field line (or rather its equiv¬ 
alence class) by its value of /i 2 ^ 2 - Furthermore, if, at a 
later time, the FLEC with this value of / 12 A 2 is found at 
a different spatial position we may define this motion as 
a field line motion knowing that this motion is uniquely 
determined by the physics because there is no remaining 
time dependent gauge freedom that changes the value of 
A2 • 

Employing the electric and magnetic fields we define the 
velocity 


B 2 + B 2 E 2& B 3 ~ hBi) 

d(h 2 A 2 ) V(h 2 A 2 ) 
dt (V(h 2 A 2 )) 2 ’ 


( 2 - 2 ) 


where ei and e 3 are unit vectors in the x\ and x 3 direc¬ 
tions, respectively, it is easy to see that w satisfies 

d( - h lf 2) + w ■ V(fe 2 A 2 ) = 0 (2-3) 

at 

and may, therefore, be regarded as the velocity, perpendic¬ 
ular to the x 2 direction, of the field lines equivalence class, 
defined by their value of h 2 A 2 . It is also clear that this 
velocity is identical to the E x B drift velocity whenever 
E-B = 0. 


2.2. Binding the particle to a field line (Equivalence 

class) 

From this point the argument is nearly the same as JKG, 
and it is straightforward to see that a particle is “bound” 
in some sense to a particular value of / 12 A 2 and hence 
to a field line. Since X 2 is an ignorable coordinate the 
X 2 component o f the canonical momentum is conserved 
( Goldstein 1980| ), h 2 P 2 + § h 2 A 2 = const. Therefore, since 
—p < P 2 < p, we have 

2cp — 

{h 2 A 2 )Max — 0^2 A 2 )Min < -^2 , (2-4) 

e 

where J 12 = [(/i 2 )Max + (^2)Min]/2. In situations where the 
particles energy and hence p ( = a/ h\p\ + h\p\ + h 2 p 2 ) 
do not change this defines a flux tube, in which the particle 
is constrained to remain. 

To obtain an estimate of the spatial scale of this flux 
tube consider a curve, S , in the x 3 ,x 3 plane, from the 
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field line with the minimum value of h 2 A 2 to the one with 
the maximum value, that is everywhere perpendicular to 
the field B'{x 1 , 13 ) = eiBi(xi, X 3 ) + e 3 B 3 (x\, x 3 ) where 
ds ■ B' = 0 along the curve S. We have then 


{h 2 A 2 ) Max (^2-^.2 )Min 



d(h 2 A 2 ) 

dx± 



d(h 2 A 2 ) 

dx 3 


(2-5) 


is clearly the azimuthal component of the canonical mo¬ 
mentum. In addition, the motion of the particles in the 
remaining degrees of freedom can often be understood by 
viewing the conserved momentum in the Hamiltonian as 
part of a potential function (Stern 1975). In such a case, 
the actual symmetry of the system aids in the understand¬ 
ing of particle dynamics. 

3.2. Problems in Charged Particle Transport 


= - [ (h 2 dsxB ') 2 > (h 2 )Min / ds \B'\ 

Js Js 

employing Equation (A2). Since the absolute value of this 
integral < — h 2 we may divide equation ( | 2 - 6 | ) by this 
quantity to obtain 




{h 2 )Max. 

(^Min 


( 2 - 6 ) 


where r' g = cp/eB' is the local gyroradius of a particle 
due to the X\ and x 3 components of the magnetic field. 
A particle, therefore, moves a distance perpendicular to 
the Xi , £3 component of the order of the harmonic mean 
of its gyroradius in this field along this path. Clearly, the 
electromagnetic fields may be arbitrary, consistent with 
Maxwell’s equations and reduced dimensionality. 

It should be noted, however, that if B lies entirely in 
the x 2 direction over a region of space, d(h 2 A 2 ) / dx 3 and 
d(h 2 A 2 )/dx\ are zero and r' g —* 00 . There is, therefore, 
no restriction, imposed by the theorem, on particle motion 
in this region. 


3. SPECIFIC CONSEQUENCES 

The theorem established in the previous section makes 
it clear that a full discussion of any problem involving 
charged particles and electromagnetic fields must either be 
done in three dimensions or the effects of reduced dimen¬ 
sionality must be shown to be irrelevant to the problem 
at hand. This has important consequences for a number 
of problems, including the physics of collisionless plasmas 
and cosmic-ray transport. We consider briefly some spe¬ 
cific applications involving both physically real and artifi¬ 
cially constrained symmetries. 


3.1. L Shell Conservation 


We will begin by considering a case in which an iden¬ 
tifiable symmetry reflects (at least approximately) a real 
property of the system. It has been known for many years 
that the particles trapped in the geomagnetic field (the 
Van Allen belts) drift around the earth in such a way as 
to approximately conserve the “L” parameter. This pa¬ 
rameter is essentially a label for the FLEC generated by 
the rotational symmetry of the dipole field. Since this 
symmetry is broken in many ways (the dipole is tilted 
with respect to the earth’s rotation axis and many per¬ 
turbations of the field, both time dependant and static, 
do not preserve this symmetry) the conservation is only 
approximate but is good enough to be a useful concept in 
space physics. Furthermore the entire concept of trapped 
(radiation belts) and excluded particles (c osmic-ray geo- 
magnetic cutoff) was discussed by Stormer ( Btormer 1955[ ) 
in terms of the “second first integral” of the motion, which 


Previous discussions of co smic-ray tran sport using the 
quasi-linear approximation ( Jokipii 1966| ) often assumed 
irregularities which depended only on one coordinate (slab 
turbulence). In such geometries the field-line mixing, or 
random walk was the only contribution to motion nor¬ 
mal to the average magnetic field. Yet, other studies of 
phenomenological particle scattering in a m agnetic field 
(Forman, Jokipii & Owens 1974, Jones, 1990) showed that 
cross field diffusion should occur. It is now clear that the 
reduced dimensionality of the slab model, not the quasi- 
linear approximation, guaranteed this result and si nce real 


turbulence is usually three dimensional in nature ( Bieber^ 
Wanner fe Matthaeus 1996| ) this limitation is probably not 
real, although the field line mixing probably plays an im¬ 
portant role. 

Recently Giacalone and Jokipii (1994) considered a 
number of issues in the general problem of cosmic- 
ray transport and its relation to the dimensionality of 
magnetic-field fluctuations. Using synthesized magnetic- 
field turbulence having a specific power spectrum, they 
directly verified that particles cannot cross field lines if the 
spectrum does not depend on all three spatial variables. 
However, when the turbulence was allowed to be fully three 
dimensional, cross field diffusion of the particles was ev¬ 
ident. This has profound implications for simulations of 
particle acceleration by collisionless shocks. 


3.3. Collisionless Shocks 

As discussed in the introduction, both hybrid and 
full particle one- and two-dimensional simulations do not 
predict any populations of ions accelerated to highly 
suprathermal energies in quasi-perpendicular shocks, This 
provides a major problem for these approaches: they are in 
direct conflict with ubiquitous observations of such popu¬ 
lations near highly oblique traveling interplanetary shocks 
(e.g. Baring et al. 1997), and the general consensus that 
the anomalous cosmic ray populations result from the ac¬ 
celeration of pick-up ions at the solar wind termination 
shock (e.g. Pesses, Jokipii & Eichler 1981). The discus¬ 
sion of JKG pointed out that the important problem of 
charged-particle acceleration at quasi-perpendicular colli¬ 
sionless shocks could not be properly studied using one- or 
two-dimensional hybrid simulations. This is because cross¬ 
field motion plays a crucial role. The restrictions for the 
acceleration at a purely perpendicular shock are particular 
severe. Here the magnetic field is carried downstream at 
the post-shock flow speed. If the particles are to remain 
tied to a magnetic field line, as required by the theorem, 
any acceleration at the shock would have to be fast enough 
for the particle gyro-radius to increase as fast as the dis¬ 
tance of the magnetic field line downstream. Otherwise 
the particle would have to be convected downstream with 
the fluid, preventing any further acceleration. It is readily 
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shown that the ordinary gradient drift along the shock face 
is too slow to effect the required increas e in the gyroradius; 
however, see the discussion in Section 3.4. Hence signifi¬ 


cant particle acceleration depends critically on the charged 
particles being able to scatter and diffuse across the mag¬ 
netic field, so that they may continue interacting with the 
shock. This necessary motion is artificially suppressed in 
simulations of reduced dimensions. Therefore one would 
incorrectly infer that significant acceleration cannot oc¬ 
cur in such quasi-perpendicular shocks. The situation in 
oblique shocks is not so severe since particles may be able 
to follow field lines to continue interacting with the shock. 
But in such cases the quantitative character of the results 
will be affected because the cross-field transport is still 
important for efficient injection. 

Furthermore, simulations of quasi-parallel shocks, where 
it has been noted that acceleration occurs by drif ting along 


a tra nsverse flu ctuation in the magnetic held (Kucharek 
& Scl joler 1991 ), should also be re-examined using three- 
dimensional simulations. In fact, Kucharek, Fujimoto & 
Scholer (1993) demonstrated that significant changes are 
incurred in quasi-parallel shock simulations in going from 
one-dimensional to two-dimensional constructions. 

3.4. Acceleration by Field Line Merging and Shock 
Surfing 

As examples of situations in which the conservation of 
the canonical momentum component does not prevent par¬ 
ticle acceleration we examine two cases where this conser¬ 
vatio n can actually cause the acceleration. There exist 
situat ions in which a charged particle is prevented from 

remaining near to the FLEC on which it was once found. 
This can arise because the FLEC ceases to exist (held line 
merging) or a strong electric held temporarily prevents 
the p article from following the held line motion (shock 

surhn b). 


For the hrst example consider a planar current sheet in 
the y — z plane which separates two regions of oppositely- 
directed magnetic held, B = e*-HoSgn(a;), where Bq is a 
constant. We can choose A y {x) = B 0 \x\ so that held lines 
that are equidistant from the neutral sheet but on oppo¬ 
site sides (therefore of opposite sign) will have the same 
value of A y . A given particle in the vicinity of this neutral 
sheet could range a further distance in x than it could if 
the held did not change sign but must still maintain the 
constancy of the quantity p y + ^A y . 

If now, as shown in Figure |l], the regions of oppositely 
directed helds are allowed to merge and annihilate held 
through plasma processes, which we do not discuss here, 
the value of A y will everywhere continuously increase due 
to the constant loss of the smallest values in the vicin¬ 
ity of x = 0. In order to preserve the invariant value of 
p y + | A y , the value of \p y \ must continuously increase. 
Thus p must continuously increase and the particle must 
be accelerated by the annihilating helds. Amb rosiano, et 
al. (Ambrosiano, Matthaeus & Goldstein 1988) employ an 
argument similar to the above in a discussion of particle 
acceleration. 
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Fig. 1.— Motion of field lines in a region of field line merging 
(annihilation) showing the increase of the local value of A y . 

In a study of plasma shocks wi th a cross shock e lectric 
potentia l drop it has been shown (Zank et al. 1996, Lee ei 


al., 199q) that some particles that are incident on the shock 


front will not have enough kinetic energy to overcome this 
potential and thus can not cross the shock because the 
sign of the potential dro p is always in the direction to re¬ 


tard the incom ing ions (Goodrich & Scudder 1984, Jonei. 
fe Ellison 1987|) . Once again, if there is an ignorable direc¬ 


tion perpendicular to the magnetic field the particle must 
remain within a gyroradius or so of it’s original FLEC, 
which is proceeding through the shock front without the 
particle. The only way this can occur is for the parti¬ 
cle’s gyroradius to increase as fast as the field lines move 
downstream. This process, called shock surfing because 
the particle skims along the shock surface along the ignor¬ 
able direction, continues until the particle gains sufficient 
energy to overcome the cross shock potential. 

Since both of these mechanisms depend rather strongly 
on there being no variation of the electric and magnetic 
fields in the ignorable direction any perturbation that does 
not have that symmetry can “derail” the process and limit 
the energy gain to one considerably less than the theoret¬ 
ical limit based on that symmetry. Therefore, once again, 
theoretical investigations of these processes should allow 
for variations in three dimensions to determine the true 
limits on their effectiveness as accelerating mechanisms. 

4. CONCLUSIONS 

In this paper we have provided a rigorous proof of the 
JKG theorem and emphasized that the most important ap¬ 
plication of the theorem is that the assumption of reduced 
dimensionality often used in modeling charged-particle 
motion has severe effects on the validity of the conclu¬ 
sions which may be drawn. In particular, any processes 
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which may depend on particle motion across the magnetic 
held cannot in principle occur. This therefore casts doubt 
on those models or simulations which are complex enough 
that the contribution of cross-held motion cannot be inde¬ 
pendently assessed. In such cases the only recourse is to 
do a fully three-dimensional model. 
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APPENDIX 


PROOF THAT A SUITABLE GAUGE EXISTS 


Since there exist many symmetries other than linear translation for which electromagnetic fields are invariant we consider 
an arbitrary, curvilinear, orthogonal coordinate system with coordinates x\, X 2 and X 3 . In general such a coordinate system 
the line element is given by 

ds = e\h\dxi + e 2 fi 2 dx 2 + e 3 h 3 dx 3 (Al) 

where the e* symbol is the unit vector in the ith direction and h\ , / 12 , and h 3 are the length scale factors for the three 
coordinates (for a Cartesian coordinate system hi = /i 2 = h 3 = 1.) In this system we have 


B = V xA = 


1 


— Zi2e2 


hih^hs 
( 9(h 3 A 3 ) 

V 9xi 


nei 


f 9(h 3 A 3 ) 9(h.2A2) 


V 9x2 


d(hiAi 
9x 3 


h 3 e 3 


9x 3 
9{h2A2) 
dx\ 


9{hiAi) \ 

9x2 ) 


(A2) 


Consider an electromagnetic field E(r,t), B(f,t) in which a particle of charge e and rest mass mo moves with velocity v. 
The only restriction is that E and B are independent of one of the spatial coordinates, say X 2 ■ 

We first demonstrate explicitly that if the electromagnetic field is not a function of a given coordinate the vector 
potential A may be chosen to be independent of this coordinate as well. The first point to note is that this independence 
is manifest by the fact that the Held structure may be moved an amount ds in the direction of this ignorable coordinate 
(if the coordinate is an angle variable the motion is a rotation) and the resulting field will be identical to the original. 
We emphasize that the transformations that follow are not coordinate transformations; they are translations or rotations 
of the fields, potentials and gauge fields themselves. The transforms are not applied to Maxwell’s Equations, rather, 
Maxwell’s Equations are applied to the transformed fields. Therefore an infinitesimal transformation gives 


where ds = \ds\. Since 


SB ' 

1 , , - 


1 ds d 



> = ds e 2 ■ V< 



1 \ 

5E\ 

l A 

| /l2 9X2 

IAI 

<L4) 


r A) 

ds d I 

'A 

H J 

) = ds e 2 ■ V < 

,9} 

h2 9X2 1 

a) 


7^0 


(in general) 


6 B = B' - B = V x (A' - A) 
= V x ( 6 A) 


(A3) 


(A4) 
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and 


we have 


A' = A + ds 


dA 
hodx 2 


- -.1 dA 

SB = ds V x —— = 0, 
h 2 ox 2 

so we may express the partial derivative with respect to x 2 as a potential field because it’s curl vanishes; 

1 dA 

-= — 

hi dxi 


where 4' is a scalar function. Since, similarly, for E we have 


1 d 


the condition on E is 


SE = E’ -E = -V(<50) - - — (SA) = 0, 

c at 


1 dE _ ^ 1 dcj) 1 d 1 dA 
hi dxi hi dxi c dt hi dxi 


= -V 


1 d(f> 1 chi' 

hi dxi c dt 


= 0 . 


This requires 


where F(t) is a function of t alone. 

We now may make a gauge transformation 


where A is an arbitrary scalar function. 
If we choose 


A = 


1 d( t> _ =F(t) 

h 2 dxi c dt 


A —► A' = A + VA 

, ,, , 1 dA 

0^0 =<P~ TT7 

c dt 


rx 2 px 2 

/ A>hi dx' 2 + K(t) / h 2 dx' 2 
Jc Jc 


(A5) 

(A6) 

(A7) 

(A8) 

(A9) 

(A10) 

(AH) 


(A12) 


(A13) 


where f^ 2 is an integral running from a fixed surface with xi = const, along a curve of constant Xi,X 3 and I\(t) is an 
arbitrary function of time alone, we see that 


1 dA 


1 dA 




1 dX 


h 2 dxi hi dxi 

and we have gauged away any xi va riatio n of A'. 


hi dxi 


= V — 


1 dX 


hi dxi 


-At = 0 


(A14) 


Similarly we have using equations All and A13 


1 dfl _ J_ d0 _ 1 d_ J__dA_ _ }_dj>_ _ 1 _ 1 dK(t) 1 dif(t) 

h 2 9 cc2 h 2 dxi c dt hi dxi h 2 dxi c dt c dt c dt 


so if we choose 


K(t) = J cF(t')dt' 


I<o 


(A15) 


(A16) 


where Kq is an arbitrary constant, we may also gauge away the x 2 variation in 0 as well. 

We can readily see that the value of A 2 is fixed, up to an overall constant in space and time, by the requirement that 
both A and 0 are independent of the Xi coordinate. For if one performs a further gauge transformation 


A' = A + VA 

0' = 0-I^ 

c dt 


(A17) 

(A18) 
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while requiring 


we immediately obtain the requirement 


t- 77 — = 7 — 77 — = 0 and 




i dA 

1 dA' 

h 2 dx 2 

h 2 dx 2 

1 d<j> 

1 d(j>' 

h 2 dx 2 

h 2 dx 2 

v( 

' 1 dA 

V l 

, h 2 dx 2 

1 d 1 

' 1 dA 

cdt \ 

, h 2 dx 2 


= 0 


= 0 


and h 2 1 dA/dx 2 (= A 2 — A 2 ) is a constant in space and time, which is merely a new choice of Kq. 


(A19) 


(A20) 



